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ABSTRACT 

There are few observational constraints on how the escape of ionizing photons from starburst galaxies 
depends on galactic parameters. Here, we report on the first major detection of an ionization cone 
in NGC 5253, a nearby starburst galaxy. This high-excitation feature is identified by mapping the 
emission-line ratios in the galaxy using [S III] A9069, [S II] A6716, and Ha narrow-band images from 
the Maryland-Magellan Tunable Filter at Las Campanas Observatory. The ionization cone appears 
optically thin, which is suggestive of the escape of ionizing photons. The cone morphology is narrow 
with an estimated solid angle covering just 3% of 47r steradians, and the young, massive clusters of the 
nuclear starburst can easily generate the radiation required to ionize the cone. Although less likely, we 
cannot rule out the possibility of an obscured AGN source. An echelle spectrum along the minor axis 
shows complex kinematics that are consistent with outflow activity. The narrow morphology of the 
ionization cone supports the scenario that an orientation bias contributes to the difficulty in detecting 
Lyman continuum emission from starbursts and Lyman break galaxies. 

Subject headings: galaxies: evolution — galaxies: individual (NGC 5253) — galaxies: ISM — galaxies: 
starburst — ISM: jets and outflows — radiative transfer 



1. INTRODUCTION 

The fate of ionizing photons in starburst galaxies is 
a critical issue for our understanding of cosmic reion- 
ization. Early-epoch star-forming galaxies are the most 
likel y candidate source for the required ionizing radiation 
(e.g., |Madau et al.|1999||Fan et al.|200T| ). However, quan- 
tify ingTFie"escapeiracHoirorTo^ (/esc) has 
proven to be challenging. 

Generally, the interstellar medium (ISM) is optically 
thick t o Lym an continuum. As discussed by IHeckman 
et al. (20111, the average hydrogen column density in 
galaxies ranges from a factor 10^ to 10^ times higher 
than the column that produces unity optical depth in 
the Lyman continuum. Thus, /osc will be dependent on 
the morphology of the ISM. For example, models have 
shown that a clumpy ISM yiel ds a higher /pa,, than a 
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smoothly v arying medium (e.g., Ciardi et al. 2002 Fer 
nandez fc Shull„2011) . For ionizing radiation to escape 
low density patns out of the galaxy must be created. In 
starburst galaxies, the mechanical feedback from massive 
stars creates low density bubbles in the ISM, which can 
then go on to break out of t he galactic disk, f acilitating 
the escape of radia tion (e.g., Mac Low & McCray 1988] 
Clarke fc Oey| |2002). However, the influence ot feedback 
on /esc is not straightforward. Simulations show that the 
very shells and bubbles that generate low-density holes in 
the ISM can initially trap ionizing ph otons be fore break- 
out, which will affect /esc (Fujita et al. 2003 ; ,Dove et al. 
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Understanding this problem is complicated by the 
small number of galaxies detected with significant /osc, 
despite large effort by the community to measure it both 
locally and at high redshift. Direct detection of excess 
Lyman continuum has bee n found in a subse t of z ^ 3 
Lyman-break galaxies (e.g. Steidel et al. 2001 1. However, 



the fraction of g alaxies with high fesc is low, on the or 



der of 10% (e.g., |Iwata eTaI]|2009| [Shapley et allp O'e 



For local starbursts, the detection rate drops, and most 
studies obtain upper limits to /esc of a few percent (e.g., 
ILeitherer et al.lll995 [ [ Heckman et al.||2001 Grimes et al. 



2009; Siana et al. 20Tof~ 

One possible explanation for the low detection rate 
is an orientation bias. In the paradigm outlined above, 
feedback processes will create passageways perpendicu- 
lar to the galaxy plane. It follows that ionizing photons 
will preferentially escape in that direction with an 
opening ang l e dependent 



2011[ ). Direct measurements of /e 



tic conditions ( 
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the orientation of the galaxy to our line-of-sight. This 
model is conducive to the formation of photoionized 
ionization cones. While photoionized emission ha s been 
observed in outflows near starburst nuclei 
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fc Bland-Hawthorn| [20101 [Shopbell fc Bland-Hawthorn 



1998[ |, clear evidence for ionization cones extending well 
into the galaxy halo has not yet been observed. Here, 
we report a new detection of such an ionization cone 
in a starburst galaxy, NGC 5253. Its narrow opening 
angle suggests that orientation may play an important 
role in the detectability of escaping Lyman continuum 
radiation. 



2. OBSERVATIONS 

NGC 5253 is a nearby, dwarf galaxy undergoing in- 
tense, centrally concentrated star formation. The most 
recent episode of star formation produced massive super 
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star clusters with ages of just a few Myr (e.g. , iCalzetti] 
et al.|1997| . At a distance of 3.8 Mpc ( |Sakai et al. 2004), 



we are able resolve structure down to tens of parsecs 
NGC 5253 is well studied across many wavelengths, pro- 
viding much ancillary data to augment this study. 

To identify channels where Lyman continuum radia- 
tion might escape from the galaxy, we apply the tech- 
nique of ionization parameter mapping (e.g., Pellegrini 
et al., in preparation; Poggc 1988). The ionization pa- 
rameter is a measure of the ionizing radiation energy den- 
sity relative to the gas density, and it can be mapped by 
emission-line ratios that probe high- vs low-ionization 
species. Regions that are optically thick to ionizing radi- 
ation show a low-ionization boundary between the highly 
ionized regions close to the ionizing source and the neu- 
tral environment; whereas optically thin regions main- 
tain high excitation throughout. For this study, we use 
[S III] A9069 and [S II] A6716 to map the high- and low- 
ionization species, respectively. 

We obtained narrow-band images of NGC 5253 on the 
nights of 2009 July 9-1 1 using the Maryla nd-Magellan 
Tunable Filter (MMTF; |\^illeux et al.|2010[ ). The band- 
passes include Ha and narrow- band continuum imaging 
observed at the redshift of NGC 5253. They are given in 
Table[TJ together with their exposure times. MMTF uses 
a low-order Fabry-Perot etalon and is mounted on the 
imaging spectrograph, IMACS, of the Magellan Baade 
telescope at Las Campanas Observatory. The IMACS 
f/2 camera has an 8K x 8K Mosaic CCD with a pixel 
scale of 0.2" per pixel. For these observations, MMTF 
was used at low etalon spacings to provide a monochro- 
matic field of view of 11.5'. We additionally obtained a 
long-slit spectrum of NGC 5253. The seeing conditions 
ranged from 0.5"-0.7" for the [S II] and [S III] images 
and 1.5"-2" for Ha. 

The bias subtraction, flat-fielding, and sky subtraction 
are accomplished using version 1 .4 of the MMTF dat a 
reduction pipeline described in Veilleux et al. (2010). 
The Ha-continuum image is first hux cahbrated using 
photometry from standard star observations. Then, the 
long slit spectrum is corrected for instrument response 
using standard star spectra, and flux calibrated using the 
flux-calibrated, Ha-continuum image. Finally, the other 
emission-line images are normalized to the emission mea- 
sured across each bandpass in the flux-calibrated, long 
slit spectrum. Systematic errors for the flux calibration 
are about a factor of two. The images are registered 
and rescaled using standard IRAfQ tasks, then contin- 
uum subtracted. Galactic extinction i s corrected assum- 
ing E{B-V) = 0.0475 across the fi eld ([Burstein fc Heiles^ 
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19821) and the extinction law from Cardeili et al 



I'he internal extinction is known to be variable t 



out the galaxy (Calzctti et al. 1997 



1989 Cresci et al. ,2005i), and its ettects are discussed in 



irough 



Caldwell & Phillips 



yx. 

A three-color composite of our data is presented in Fig- 
ure [l] along with the continuum-subtracted, emission- 
line images. Sky background and noise properties of the 
images are measured using the median and standard de- 



® http:/ /www. astro. umd.edu/'veilleux/mmtf/datared. html 
IRAF is distributed by NOAO, which is operated by AURA, 
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dation. 



viation in the flux of fifteen 20" x 20" boxes around the 
galaxy. In Ha, we have a 2>a confidence detection limit of 
4.1 X 10""'^^ ergs""'^ cm~^. The ionized gas in the galaxy 
has a roughly spherical distribution around the nucleus, 
with networks of loops and filaments par ticularly notice- 
able to the south and northwest (e.g., Marlowe et al. 
19951). 






Figure 1. a) Composite image of NGC 5253. Red, blue and green 
correspond to [S III], [S II], and continuum at A6680, respectively. 
Ionization cone extends ESE from the nuclear starburst. Panels h, 
c, and d show individual continuum-subtracted, [S III], [S II], and 
Ha, respectively. At a distance of 3.8 Mpc, 10"= 180 pc. 



3. IONIZATION CONE IN NGC 5253 

3.1. Morphology and Excitation of the Ionization Cone 

NGC 5253 shows a stunning detection of an ionization 
cone in a starburst galaxy. In Figure [T] the ionization 
cone stands out beautifully, extending ESE along the mi- 
nor axis. In the ratio map. Figure [2] the ionization cone 
is identified by its high excitation relative to the rest of 
the galaxy. This high excitatio n has also been observed 
in [O III]/[S II] dCalzetti et aI][T999i ). It has an opening 
angle < 40° with a narrow appearance at large radii. As 
the distance from the galaxy increases, the morphology 
of the cone becomes filamentary. While the cone appears 
to be one-sided, we note that extended [S HI] emission 
to the northwest may suggest a bi-cone (Figure [T]). 
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Figure 2. [S III]/[S II] ratio map for NGC 5253. Dark indicates 
high values of [S III]/[S II]. 

The changes in ionization parameter as a function of 
radius are shown more quantitatively in Figure |3] These 
plots are generated using 3x3 binning of the continuum- 
subtracted images, and a 3a detection cut-off in all three 
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Table 1 

Observations 



Emission 


Wavelength 


Effective 


Exp. Time 


Flux=» 


line 


(A) 


Bandpass (A) 




(10^^"^ erg s^^ cm^^ ) 


[S III] 


9081 


27.0 


5 X 1200 s 


0.39 


continuum 


9180 


27.0 


5 X 1200 s 




[S II] 


6726 


16.2 


5 X 1200 s 


0.59 


continuum 


6680 


16.2 


5 X 1200 s 




Ha 


6572 


18.5 


3 X 1200 s 


12.9 


continuum 


6680 


18.5 


3 X 1200 s 




^ Continuum- 


subtracted flux measured in 


tlie ionization cone using the DS9 funtools module. 
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Figure 3. Emission-line ratios vs. galactocentric radius. Red 
represents ionization cone emission. Blue line is the mean galac- 
tic value in radial bins, vfith bars noting one standard deviation. 
Representative measurement error is shown by the green error bar. 
The black lines indicate the effect of different levels of extinc- 
tion {Ay = 1,2,5,10). Top, middle, and bottom panels show 
[S III]/[S II], [S II] /Ha, and [S HI] /Ha, respectively. 



emission lines. The ionization cone, shown in red, is 
defined as the region between P.A. — 102°-140° using 
a = 13'^39'"55^7, (5 = -31°38'24" (J2000) as the apex. 
This point, marked by the cross in Figure [2j is the lo- 
cation of a likely source for the ionizing radiation, as 
discussed below. For comparison, the blue line is repre- 
sentative of the rest of the galaxy. 

The top panel of Figure ^ clearly shows a strong ex- 
cess in the [S ni]/[S II] ratio for the ionization cone at 
lower galactocentric radii R, and a clear excess that con- 
tinues to i? ~ 600 — 700 pc. The middle and lower panels 
show that this excess cannot be explained solely by the 
well-known increasing gradient in [S II] /Ha that is found 
in the diffu se, ionized gas of star-forming galaxies (e.g^ 
Rand||1998[ [Haffner et al]|1999| ). Furthermore, Figurejl] 
shows that the primary star-forming region extends only 
to ~ 500 pc, whereas photoionized material associated 
with the cone clearly extends beyond. This strongly sug- 
gests that there is not enough gas in the outer halo to 
absorb all the ionizing photons, and that they therefore 
may escape the galaxy. 

As noted abov e, the extinction t hroughout the star- 
burst is variable ( Cresci et aH2005 ), and might enhance 
the observed [S IiiJ/[y iij ratio relative to the intrinsic 
one. In particular, a dust lane with associated CO along 
the minor axis is coincident with the location of the cone 



al. 



2002|. However, Figure 3 shows that we need 



et 

Ay ~ 3— lO mag to explain the observed [S III]/[S II] ra- 
tios entirely by the dust lane, which is m uch more than 
its ob served extinction {Ay ^ 2.2 mag, Calzetti et al. 
If there were that much extinction, we would ex- 
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pect a much lower flux along the cone, which is not seen 
in the emission-line images in Figure [T] Additionally, the 
[S II] /Ha ratio observed in the cone cannot be explained 
by extinction. We also note that the highest extinction 
will likely be strongest at the smallest R. 

Another possible effect is a bias in our observed 
[S in]/[S II] ratio caused by strong variations in electron 
density. This affects the [S II]A6716/A6731 ratio, so that 
our observed A6716 may underestimate the total [S II] at 
high density. This error again will be largest at smaller 
radii where densities are larger, and leads to difference 
of ~ 25% between a density of 10^ and lO'^ cm~'^. In our 
long-slit spectrum, [S II] AA6716/6731 ratio ranges from 
1.4 to 1.0, which correspond to densities of < 100 cm~'^ 
and ^ 500 cm~^, respectively. Since this is a smaller 
variation than between 10^ and 10^ cm""^, the error on 
our observed [S II] due to density variations is likely much 
less than 25%. 

The measured emission-line fluxes in the ionization 
cone (Table [l]) can set limits on the ionizing source. 
The Ha flux corresponds to an Ha luminosity L(Ha) — 
2.2(±l.l)x lO'^Q ergs-^ This is ^ 6% the total L(Ha) in 
the galaxy ( Marlowe et al.|[l995 corrected for distance), 
from which we obtain a lower limit of log(Qo) > 51.19 
on the rate of ionizing photons, Qo{s~^), needed to 
explain the observed emission. A n 07 star with the 
- 0.2 ^0 metalhcit y of NGC 5253 (|Walsh fc Roy||1989 l 
has log(Qo) = 49.0 ( [Smith et al.|2002l ). I'his imphes that 
radiation equivalent to that produced by > 160 ± 80 07 
stars is needed to ionize the gas. This is a lower limit for 
the radiation escaping the nuclear starburst, since the 
flux will be affected by internal extinction and /osc- 

NGC 5253 is host to many young clusters that could 
be the source for the ionizing radiation illumin ating the 
cone. Two likely candidates are NGC 5253-5 (Calzetti 
and the ionizing cluster of the radio su- 
Turner et al.'''2000; Gorjian 



et al. 



perne 



et al. 



1997 



3ula dBeck et al. 



1996 



2001 ). I'he location o: 



these objects, separated by 



(Calzetti et al. 1997 Caldwell & Phillips 1989 Meier 



just a few arcseconds on the sky ( Alonso-Herrero et al. 
20041, is noted by the cross in Figure 2 Both clusters 
are very young and very massive d Turner et al. 2000 
Calzetti et al.|1997[ [Alonso-Herrero et af.|2U04[ ). In fact, 
the radio supernebu fa may be gravitationalfy confined 
(Turner et al. 2003). The 20 pc region around the ra- 
dio supernebufa hosts an ionizing population ~ 7000 07 
stars, with some uncertainty due to the high extinction 
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Figure 4. NGC 5253 Ho image, minor axis oriented along the 
y-axis. Ha echelle profiles are overlaid to the left. The profiles are 
normalized to a peak of one and the spacing between each profile is 
2. "6. The red l ine denotes the Vsys fo r NGC 5253, reported to be 
407ifc3 kms-l |Koribalski et al.|2004[ . The slit position is marked 
in yellow. 



towards this cluster (Turner & Beck 2004). Thus, the 
clusters in NGC 5253 easily generate enougli ionizing flux 
to explain the ionization cone. 

3.2. Kinematics 

The HI kinematics in NGC 52 53 show only small scale 
rotat ion about the minor axis (Kobulnicky & Skillman 
2008| ), making it unh kely that NGC 5253 is rotation- 
allysupported ( jCaldwell fc Phillips, 1989). In Fig ure 
the line profiles trom an Ha echellegram (Slit 3, Mar- 
tin & Kennicutt 1995) are overlaid on the H« image ot 
NGC 5253. 'i'he velocity shows significant variation with 
position, and a maximum projected velocity difference 
Av ~ 30 kms^^ relative to systemic. 

An expanding wind might explain the kinematics. In 
this scenario, we expect the velocity to increase with dis- 
tance from the galaxy, since the wind is expanding into 
the lower density gas of the halo. To the southeast, we do 
see an increase in velocity with galactocentric distance. 
However, it reaches a maximum at ^ 40" from the galaxy 
center before decreasing towards systemic at larger R. 
This does not strictly follow what is expected for a free- 
flowing wind, in which the velocity should continue to 
increase. Additionally, an expanding wind typically ex- 
hibits line splitting, which is n ot seen in the Ha pro- 
files (Martin & Kennicutt 1995). Kinematic data from 



HI observations are sugg estive of either inflow or out - 
flow along the minor axis ( [Kobulnicky fc Skillman(2008| . 
Therefore, although not completely straightforward, me- 
chanical feedback might plausibly explain the observed 
kinematics based on the limited information available. 



3.3. Possible AGN 

Typically, ionization cones are associated with AGN, 
which suggests an alternative explanation for the ioniza- 
tion cone in NGC 5253. The morphology of the cone 
is very narrow. This suggests coUimation of the ioniz- 
ing radiation, as would be expected if there is an AGN 
present. The radio supernebula is particularly intrigu- 
ing in this context. Only 0."6 in size, it provides 
roughly 80% of the galaxy's IR flux at 12^m (Gorjian 



et al. 2001). There is some debate about whether a 



young super star cluster or AGN is the ionizing source 
of this nebula. Its radio contin uum has a flat slope, con- 
sistent with thermal emission (Beck et al. 1996 Turner 



et al.||1998[ |2000 [), but which do es not rule out an AGN 
component ( I'urner et al.||1998 |. However, the observed 



spectral lines ar e also narrower th an would be expected 
from an AGN ( [Beck et ar]|1996[ ). Additionally, X-ray 
observations of NGC 5253 show that the point sources 
and diffuse emission are consistent with star formation 
and shocks (Ott e t al.l2005[ [Summers et al.|2004[ [Martin 
[fc^ Kennicutt 1 9U5 ). 'iheretore, while we cannot rule out 
an AGN m NgC 5253, there are challenges facing that 
interpretation. 

If NGC 5253 has an obscured AGN, a precessing jet 
might explain the gas kinematics. In this case, the ve- 
locities on one side of the galaxy should be mirrored ex- 
actly on the opposite side. In NGC 5253, the overall 
sinusoid in velocity to the southeast is roughly reflected 
to the northwest. Along the ionization cone, the velocity 
transitions from redshift to blueshift around ~ 30" from 
the galaxy. To the northwest, the opposite transition, of 
roughly equal amplitude, is seen at ~ 28". However, the 
reflection is not perfect; there are points to the north- 
west where the velocity briefly flips, which are not seen 
to the southeast. Additionally, the amount of blueshift 
to the northwest does not match the redshift observed in 
the southeast. Therefore, while the general kinematics 
could be explained by a precessing jet, the discrepancies 
are difficult to reconcile in detail. 

4. IMPLICATIONS: ORIENTATION BIAS 

Based on our detection of an ionization cone in 
NGC 5253, we suggest that an orientation bias can ex- 
plain the difficulty in detecting escaping ionizing radia- 
tion in starburst galaxies. We observe evidence that ion- 
izing radiation from NGC 5253 is traveling to large radii, 
if not escaping the galaxy entirely. Assuming the cone 
is axisymmetric, the estimated solid angle subtended by 
the cone is ~ 3% of 47r steradians, suggesting that the 
escape process happens over a small solid angle. If we 
assume isotropic radiation from 7000 07 stars, as above, 
we could expect to see ionized-gas emission in this solid 
angle equivalent to that produced by 210 07 stars. This 
is similar to our estimated 160 ± 80 07 stars needed to 
generate the observed Ho emission in the cone, support- 
ing the physical link between the source and the cone. 
The narrow morphology suggests that in order to detect 
Lyman continuum, the line-of-sight to the galaxy must 
be close to the axis of escape. Thus, the orientation of 
the galaxy should strongly influence the detectability of 
Lyman continuum radiation. 

We can find further evidence for this scenario in recent 



observations by Heckman et al. ( 2011, ) . In their study of 
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11 Lyman-break analog galaxies and 15 local starbursts, 
including NGC 5253, they found indirect evidence for sig- 
nificant /esc in three of the Lyman-break analogs. Out- 
flow velocities on the order of 10"^ kms"^ were found in 
all the galaxies suggested to have significant f^sc- These 
outflow velocities are so high, that we can interpret it to 
mean the direction of the outflow must be closely aligned 
to the lines-of-sight towards the galaxies. Since winds in 
starburst galaxies are generally coincident with the poles 



(Veilleux et al. 
high /esc are a! 



20051), this suggests that these inferred 
associated with face-on systems. This 



fits our explanation that measurements of large /esc are 
biased towards galaxies whose outfiow happens along our 
line-of-sight. 

5. CONCLUSION 

The ionization cone in NGC 5253 provides a new per- 
spective to our understanding of the fate of ionizing 
photons in starburst galaxies. This ionization cone ap- 
pears optically thin, which is suggestive of the escape of 
ionizing radiation along the cone. At minimum, ioniz- 
ing radiation is escaping from the nuclear starburst into 
the galaxy's halo. We considered both stellar and non- 
stellar sources for the ionization cone. Assuming a stel- 
lar source, the massive clusters in NGC 5253 produce 
enough ionizing radiation to explain the Ha emission ob- 
served. However, the possibility of an obscured AGN has 
not been ruled out. The Hq kinematics along the cone 
exhibit a complex morphology, but are consistent with 
some form of outflow. The data are not straightforward 
to interpret, but might be explained by an expanding 
wind or perhaps a precessing jet. 

Using the ionization cone of NGC 5253 as an analog 
for other starburst galaxies, we see that ionizing radia- 
tion escapes along the minor axis. Additionally, the solid 
angle over which this radiation escapes is small, based on 
the cone morphology. This will constrain the viewing an- 
gles from which Lyman continuum radiation can be de- 
tected. Thus, the orientation of the galaxy will strongly 
influence the ability to detect Lyman continuum radia- 
tion. NGC 5253 may be the nearest starburst galaxy 
with a significant /esc, and offers a unique opportunity 
to study the emission mechanisms in detail. 

We kindly thank David Osip for recovering important 
data from the Magellan Archives, John Glaspey for con- 
verting old image formats to current ones, and Eric Pel- 
legrini for useful discussions. We thank the referee for 
thorough and helpful comments. This work was funded 
by NSF grants AST-0806476 to MSG and AST-1009583 
to SV. 
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